Switching On Magnetism in Ni-doped Graphene 
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Magnetic properties of graphenic carbon nanostructures, relevant for future spintronic applica- 
tions, depend crucially on doping and on the presence of defects. In this paper we study the mag- 
netism of the recently detected substitutional Ni (Nisut) impurities, fiisub defects are non-magnetic 
in flat graphene and develop a non-zero magnetic moment only in metallic nanotubes. This surpris- 
ing behavior stems from the peculiar curvature dependence of the electronic structure of Ni^ub. A 
similar magnetic/non-magnetic transition of f^isut can be expected by applying anisotropic strain 
to a flat graphene layer. 

PACS numbers: 75.75.-(-a, 73.22.-f, 73.20.Hb, 81.05. Uw 



I. INTRODUCTION 

Graphenic carbon nanostructures have opened new re- 
search routes in nanoelectronics.-i*^ In particular, their 
magnetic properties are receiving much attention, both 
experimental and theoretically, mainly in relation to spin- 
tronics. Spin-qubits and other spintronic devices seem 
feasible due to the very long spin relaxation and deco- 
herence times in graphene^-:- and to the fact that the 
magnetism of the edge states of graphenic nanoribbons 
can be controlled by applying external electric fields i^i^ 
Magneto-optical properties are also being actively stud- 
iedi^ii All these properties are drastically influenced by 
the presence of defects and dopantsi^ For example, 
strong magnetic signals in nanocarbons have been re- 
ported after irradiation that seem to be associated with 
the creation of defects j^ ^'^^i^^i^^ Thus, the magnetic and 
transport properties of carbon systems can in principle 
be engineered using these additional degrees of freedom. 

In a recent x-ray adsorption study, Ushiro et aZ.^^ have 
demonstrated the presence of important amounts of sub- 
stitutional Ni (Nisub) impurites in purified carbon nanos- 
tructures synthesized using Ni containing catalyst. Such 
Nisub impurities were also observed by Banhart et alr^ in 
electron microscopy images of onion-like graphenic parti- 
cles. In general, due to the stability of the substitutional 
configuration, the incorporation of transition metals to 
the carbon layer during growth or saturating existing va- 
cancies seems a likely process. In spite of this, the mag- 
netic properties of substitutional transition-metal impu- 
rities in graphenic systems have not been studied in de- 
tail. Instead, most theoretical studies have focused on 
adsorbed transition- metal atoms on fuUerenes and single- 
walled carbon nanotubes (SWCNTs)^iilii^ii^i^i2i Few 
calculations to date have considered Nigub impuritie a^^i^^ 
and little attention was paid to the magnetic properties. 

In this paper we present a first-principles density func- 



tional theory (DFT) study of Nigub defects in graphene 
and armchair and zigzag SWCNTs of different diame- 
ters. We have discovered that the magnetic moment of 
substitutionally Ni-doped graphene can be controlled by 
applying mechanical deformations that break the hexago- 
nal symmetry of the layer, like curvature does. Surprins- 
ingly, Nistif, impurities are non-magnetic in flat graphene. 
However, their magnetic moment can be switched on by 
applying curvature to the structure. Furthermore, the 
magnetic moment of Nisuf, also becomes a signature of 
the metallicity of the structure: only metallic tubes de- 
velop a moment that depends on the tube diameter and 
Ni concentration. This behavior stems from the peculiar 
curvature dependence of the electronic structure of the 
Nisub impurity. 



II. METHODOLOGY 

Our calculations have been perfomed with the SIESTA 
code^i^i^^ using the generalized gradient approximation 
(GGA)^ to DFT and TrouUier-Martins^l pseudopoten- 
tials. Our Ni pseudopotential includes a pseudo-core with 
a core matching radius of 0.53 a.u. in order to perform 
non-linear core corrections for the description of exchange 
and correlationi^ We have tested that this pseudopo- 
tential yields the correct magnetic moment and band 
structure for bulk Ni. We have used an energy shift— 
of 50 meV to define the radii of the different orbitals. A 
double-C polarized (DZP) basis set has been used for the 
calculation of the magnetic moments and electronic band 
structures of all our systems. However, we have checked 
that a double-C (DZ) basis yields to almost identical re- 
laxed structures as the DZP basis and, therefore, we have 
used the smaller DZ basis for the relaxations of systems 
containing more than ~100 atoms (our supercells contain 
up to 512 atoms, see Fig. [S] below). Atomic coordinates 
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FIG. 1: (Color online) (a) Relaxed geometry of a substitu- 
tional Ni (Nisui,) impurity in a (5,5) SWCNT, and (b) iso- 
surface (±0.004 e~/Bohr^) of the magnetization density with 
light (gray) and dark (blue) surfaces corresponding, respec- 
tively, to majority and minority spin. 



were relaxed using a conjugated gradient algorithm until 
all force components were smaller than 0.05 eV/A. Rele- 
vant lattice parameters were fixed to those of the pristine 
systems. To prevent spurious interactions the minimum 
distance between the walls of neighboring SWCNTs was 
18 A. The fineness of the real-space grid used to calculate 
the Hartree and exchange-correlation contribution to the 
total energy and Hamiltonian was equivalent to a 180 Ry 
plane-wave cutoff. A good integration over the Brillouin 
zone proved to be instrumental to converge the magnetic 
moment of Nis„b impurities in metallic SWCNTs. The k- 
point sampling was equivalent to a 136x136 sampling^^ 
of the Brillouin zone of graphene. We used a Fermi- 
Dirac distribution with KbT=21 meV. In some cases 
we have checked that our results are reproduced using 
a different methodology: the VASP code with projected- 
augmented-wave potentials and a plane-wave cut-off en- 
ergy of 400 eV^^ 



III. RESULTS AND DISCUSSION 

Figure [1] (a) illustrates the typical equilibrium struc- 
ture of Nisub in the case of a (5,5) SWCNT. The Ni atom 
appears displaced ~0.9 A from the carbon plane. Al- 
though both outward and inward displacements can be 
stabilized, the outward configuration is always more sta- 
ble. The calculated Ni-C distances {d-Nic) are in the 
range 1.77-1.85 A in agreement with experimen t 
Armchair tubes exhibit two slightly shorter and one 
larger values of disnc, the opposite happens for (n,0) 
tubes, whereas for graphene we obtain a threefold sym- 
metric structure with A. Ni adsorption in- 
hibits the reconstruction''^ of the carbon vacancy. Fur- 




FIG. 2: (Color online) (a) Schematic representation of the 
electronic structure of a Nisub impurity in graphene. Panel 
(b) shows the calculated band structure for a Nisut. impu- 
rity in a 4x4 graphene supercell. The size of the circles and 
squares corresponds to the amount of Ni 3d and 4s character 
respectively. Energies are referred to the Fermi energy. 



thermore, we have checked that a symmetric structure 
is obtained even when starting from a relaxed vacancy. 
The Ni binding energy is quite large: 7.9 eV for graphene 
and about 8.5 eV for (5,5) and (8,0) tubes. The calcu- 
lated adsorption energy of Ni on the surface of the same 
tubes is ~ 2.5 eV<^ Thus, we can conclude that the for- 
mation of Nisub defects by passivation of existing carbon 
vacancies is a very likely process both for graphene and 
for SWCNTs. 

Figure [1] (b) shows the magnetization density for a 
Nisub defect in a (5,5) metallic nanotube at large dilution 
(0.3 % Ni concentration). The total magnetic moment of 
this system is 0.5 fiB- The magnetization comprises the 
Ni atom and its C neighbors. However, it also extents 
considerably along the tube, particularly in the direction 
perpendicular to the tube axis. This indicates the polar- 
ization of some of the delocalized electronic states in the 
nanotube. Indeed, as we clarify below, the magnetism in 
substitutionally Ni-doped SWCNTs only appears associ- 
ated with the curvature and the metallicity of the host 
structure 

Figure [2] (a) shows a scheme of the electronic struc- 
ture of Nisub in graphene, while Fig. O (b) presents the 
calculated band structure using a 4x4 supercell. Very 
similar results are obtained using larger supercellsj22. Sev- 
eral levels with Ni-C bonding character and a strong Ni 
3d contribution can be found between 2 and 6 eV below 
the Fermi energy (E^?). This considerable band width is 
a signature of the strong Ni-C interaction. As a conse- 
quence of the bonding interaction the Ni 3d band is sta- 
bilized and can be pictured as "fully occupied" . Close to 
Ep we find three levels with Ni-C antibonding character. 
One of them is occupied and appears around 0.7 eV below 
Ep close to r. This level comes from a fully symmetric 
linear combination of the 2pz orbitals (z-axis normal to 
the layer) of the nearest C neighbors interacting with the 
3^22 orbital of Ni. Two levels coming from the hybridiza- 
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FIG. 3; (Color online) Band structure of a (5,5) nanotube 
containing a Ni impurity every 4 unit cells (Nisut-Nisut dis- 
tance of ~9.8 A) for (a) a paramagnetic calculation, and for 
(b) majority and (c) minority spins. Circles and triangles 
correspond respectively to the amount of Ni 3dyz and 3dxz 
character. X-axis is parallel to the tube axis and y-axis is 
tangential. 



tion of the in-plane sp lobes of the carbon neighbors with 
the Ni 3dxz and 3dyz orbitals appear ^^0.5 eV above Ep. 

As a consequence of this electronic structure, with the 
Ni 3d states well below Ep and no flat bands crossing E p , 
the magnetic moment of the Ni^uf, impurity in graphene 
is zero. Interestingly, the three levels appearing close to 
Ei? in Fig. [2] are reminiscent of those found for the unre- 
constructed carbon vacancy in graphene^^. The presence 
of a Ni atom stabilizes such symmetric structure of the 
vacancy and slightly shifts the two levels coming from the 
sp lobes (now C 2sp-Ni 3d antibonding levels) to higher 
energies. This contributes to stabilize the paramagnetic 
solution and renders a Ni<,„b impurity with zero spin po- 
larization in the flat layer of graphene. 

The basic picture described above is still valid for the 
electronic structure of the Nigub impurity in SWCNTs. 
However, the modifications that appear due to the cur- 
vature of the carbon layer are responsible for the ap- 
pearance of a magnetic moment. Figured] (a) shows the 
band structure of a paramagnetic calculation of a (5,5) 
SWCNT with a Nisub impurity every four unit cells. Here 
the distance between neighboring Nisub impurities is sim- 
ilar to that of the graphene layer in Fig. [2] (b), although 
the Ni concentration is 2.5 times lower (1.3 %). Com- 
paring these two Figures we can appreciate the effects of 
curvature. The degeneracy between dxz and dyz states 
is removed (x-axis taken along the tube axis and y-axis 
along the tangential direction at the Ni site). The dyz 
contribution is stabilized by several tenths of eV and a 
quite flat band with strong dyz character is found pinned 
a.t Fip close to the Brillouin-zone boundary. Under these 
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FIG. 4: (Color online) Effect of curvature (anisotropic strain) 
on Nisub in (n,n) tubes. The electronic structure of 'Nisub is 
similar to that of the unreconstructed carbon vacancy. One of 
the impurity levels with antibonding C 2sp-Ni 3d character is 
shifted downwards and, for large enough curvatures, becomes 
partially populated and spin-polarized. 



conditions the spin-compensated solution becomes unsta- 
ble and a magnetic moment of 0.48 fiB is developed. Fig- 
ures [3] (b) and (c) show, respectively, the band structure 
for majority and minority spins. The exchange splitting 
of the dyz level is ~0.4 eV and the energy gain with re- 
spect to the paramagnetic solution is 32 meV. Similar 
results are obtained using the VASP code. 

In general, whenever a flat impurity with appreciable 
Ni 3d character becomes partially filled we can expect 
the appearance of a magnetic moment. The population of 
such an impurity level occurs at the expense of the simul- 
taneous depopulation of some of the delocalized carbon 
Pz levels of the host structure. For this reason the de- 
velopment of a magnetic moment is more likely for Nisub 
impurities in metallic structures like the armchair tubes. 
The crucial role of the host states also explains the delo- 
calized character of the magnetization density depicted 
in Fig. [1] (b). However, it is important to stress that 
the driving force for the formation of a magnetic mo- 
ment associated with the Nis„b impurity in a SWCNT 
is the local curvature of the carbon layer that shifts the 
energy position of one of the impurity levels downwards 
until it crosses Ep. A schematic representation of this 
phenomenon can be found in Fig. 3] where we also em- 
phasize the similarities between the levels of the Nigub 
defect and those of the unreconstructed carbon vacancy. 
At large tube diameters we must recover the limit of flat 
graphene with zero magnetic moment. 

For semiconducting tubes the situation is somewhat 
different. The d^z and dyz derived levels will remain 
unoccupied unless their energies are shifted by a larger 
amount that pushes one of them below the top of the va- 
lence band. Therefore, if the tube has a large enough gap 
the magnetic moment will be zero irrespective of the tube 
diameter. We have explicitly checked that a zero mag- 
netic moment is obtained for (8,0) and (10,0) semicon- 
ducting tubes for Ni concentrations ranging from 1.5% 
to 0.5%. The different magnetic behavior of Nisub im- 
purities depending on the metallic and semiconducting 
character of the host structure provides a route to exper- 
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FIG. 5: (Color online) Magnetic moment per Ni impurity 
for different (n,n) tubes as a function of the distance between 
periodic images of the impurity, i.e., the length of the super- 
cell. The vertical lines and the numbers at the top indicate 
the number of unit cells in a supercell of a given length. 

imentally identify metallic armchair tubes. 

Figure [5] displays the magnetic moment per Nisut atom 
for Ni-doped armchair tubes of different diameters. All 
of them present a magnetic moment that oscillates as a 
function of the tube diameter and the size of the super- 
cell used in the calculation, i.e., the Nis„f,-Nistib distance. 
For (4,4), (5,5) and (6,6) tubes the first supercell show- 
ing a non-zero magnetic moment contains four unit cells. 
For (7,7) and (8,8) tubes this minimum length increases 
up to eight and six unit cells, respectively. The appear- 
ance of a complex oscillatory pattern as a function of the 
Nisub-Nisuf) distance is easily understood if we recall that 
the magnetic moment critically depends on the energy 
position of a particular impurity level nearby Fip and 
the strong hybridization of this level with the delocalized 
states of the nanotube. Another consequence of this hy- 
bridization is the long range of the interaction between 
Nisub impurities: for example, the magnetic moment in 
(5,5) and (4,4) tubes still presents strong oscillations in a 
range of Nisub-Ni^ub distances between 20 and 40 A. Un- 
fortunately, a meaningful exploration of larger distances 
between impurities requires a methodology different from 
the ab initio supercell approach used here. 

Finally, we have examined the magnetic coupling be- 
tween Nisub impurities in the (5,5) tube. We have dou- 
bled some of the simulation cells considered above, so 
they contain two Ni atoms, and calculated ferromagnetic 
(FM) and antiferromagnetic (AFM) arrangements. Only 
when the two Ni<;„f, impurities shared a common C neigh- 
bor the AFM arrangement was favored by a few meV. FM 
configurations were the most stable in all other cases. Al- 
though more work is necessary to accurately determine 
the size and distance dependence of the effective exchange 
interaction (probably using model calculations similar to 



those in Ref. [33 ) , our results indicate that Nis„b impuri- 
ties in small diameter metallic nanotubes, where they de- 
velop a magnetic moment, can exhibit relatively large FM 
couplings that slowly decay with distance. For example, 
for Niguf, impurities at distances of 10 A in a (5,5) tube 
we estimate a Je/ / of ~9 meV. This FM interactions can 
have important implications for the experimental detec- 
tion of the curvature dependent magnetism of the Ni^ub 
impurities described in this work and its influence in the 
observation of magnetism in carbon nanotube samples. 



IV. CONCLUSIONS 

The calculations presented here show that substi- 
tutional Ni impurities in graphenic carbon structures 
present a strong covalent interaction with their carbon 
neighbors. This interaction stabilizes the 3d levels of 
Ni that appear as a completely flUed shell. As a conse- 
quence, Nisub is non-magnetic in graphene. Two unoccu- 
pied levels of the Nis„b impurity appear close toEp corre- 
sponding to antibonding combinations of Ni 3d orbitals 
and C 2sp lobes. Curvature breaks the degeneracy of 
these two levels and, under appropriate conditions, shifts 
one of them to lower energies such that becomes partially 
occupied in metallic armchair tubes. In this situation the 
system develops a magnetic moment that exhibits a com- 
plex behavior as a function of the tube diameter and Ni 
concentration. However, for semiconducting zigzag tubes 
the Nisub impurities remain non-magnetic. These results 
should be taken into account when studying the magnetic 
properties of carbon nanostructures. 

In summary, curvature can be used to drive substi- 
tutional Ni impurities between a magnetic and a non- 
magnetic state in metallic graphenic carbon nanostruc- 
tures. In particular, the magnetic properties of substi- 
tutionally Ni-doped graphene can be tuned by control- 
ling the curvature of the layer around the Nis„b defects. 
A possible way to do this is to deposit the graphene 
layer on a substrate with a small lattice parameter mis- 
match. This can cause the undulation and local bending 
of the graphene layer— i^. Alternatively, a similar non- 
magnetic/magnetic transition of Ni^ub might be obtained 
applying anisotropic strain to a flat graphene layer. In 
this regard, we are currently exploring if the formation 
of a magnetic moment can also be controlled by applying 
uniaxial strain to a flat Nis„b-doped graphene layer. 
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